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The exact dynamics of a Jaynes-Cummings model for a qubit interacting with a continuous dis- 
tribution of bosons, characterized by a special form of the spectral density, is evaluated analytically. 
The special reservoir is designed to induce anomalous decoherence, resulting in an inverse power 
law relaxation, of power 3/2, over an evaluated long time scale. If compared to the exponential-like 
relaxation obtained from the original Jaynes-Cummings model for Lorentzian-type spectral density 
functions, decoherence is strongly suppressed. The special reservoir exhibits a band edge frequency 
coinciding with the qubit transition frequency. Known theoretical models of photonic band gap 
media suitable for the realization of the designed reservoir are proposed. 
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I. INTRODUCTION 

The dynamics of a two-level system (TLS) coupled 
to an environment is of great interest [2 [2], given the 
many applications in quantum optics, quantum informa- 
tion, atom-cavity interactions, molecular dynamics, spec- 
troscopy, and solid state physics. Within the theory of 
open quantum systems [3] the environment is often rep- 
resented as a reservoir of bosons [J] , characterized by the 
corresponding spectral density function. Several stud- 
ies have been performed to investigate the dependance 
of decoherence on the explicit form of the spectral den- 
sity function. According to the system of relevance, sub- 
ohmic, ohmic, super-ohmic and Lorentzian forms of the 
spectral density have been investigated 0] , by adopt- 
ing the celebrated Jaynes-Cumming model [5j. 

Interesting results emerge from the technique of the 
resolvent operator [5J applied to a two-level atom (TLA) 
interacting, in rotating wave approximation, with a 
Lorentzian |7] distribution of field modes. The cou- 
pling constants are assumed to vary slowly for frequency 
changes. In this way, oscillating behaviors enveloped in 
exponential decays emerge in the exact dynamics. For a 
detailed report we refer to [5]. 

Another model mimicking the spontaneous decay of a 
TLA in a structured reservoir, has been introduced by 
Garraway [pj HO] and solved exactly for Lorentzian type 
distributions and a special non-Lorentzian one with two 
poles in the lower half plane. The assumption that the 
frequency range is (— oo, +00) makes it possible to char- 
acterize the exact dynamics through the poles of the spec- 
tral density in the lower half plane and the corresponding 
pseudomodes [IT]. The resulting dynamics is described 
by oscillations enveloped in exponential relaxations. 

Recently, the model adopted in Ref. [9] has been re- 
sumed to compare the exact master equation in time- 
convolutionless form with the Nakajima-Zwanzig master 
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equation, where the perturbation expansion of the corre- 
sponding memory kernel is performed |12j . 

In line with the attempt to delay the destructive effects 
of the environment on the qubit, i.e., decoherence, we 
consider the model adopted in Ref. [5] and study a reser- 
voir that is piecewise similar to those usually adopted, 
e.g., sub-ohmic and Lorentzian ones. We show the rele- 
vant changes in the time evolution induced by the spe- 
cially designed reservoir and how the decoherence process 
results to be strongly suppressed. We notice that the an- 
alytical calculations concerning the exact dynamics, are 
performed for a positive range of modes frequencies. 



II. THE MODEL 

The model considered was originally introduced by 
Garraway |9] in order to describe the exact dynamics 
of an atom interacting with a cavity field, approximated 
by a continuous distribution of modes, described by a 
Lorentzian, a linear combination of two Lorentzian and 
a non-Lorentzian distribution. Briefly, the Hamiltonian 
of the whole system reads H = Hg + He + Hi, where the 
system Hamiltonian Hs, the bath Hamiltonian He and 
the interaction Hamiltonian are given by the following 
forms: 



H s = lu o- + er_ 
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being h = 1 in the system of units adopted. In the usual 
notation a + and <r_ are the rising and lowering oper- 
ators, respectively, acting on the Hilbert space of the 
qubit, while a\ and aj~ are the creation and annihilation 
operators, respectively, acting on the Hilbert space of 
the fc-th boson, satisfying the boson commutation rule, 

dk,a k , = Sk,k', for every non- vanishing natural value 



2 



of k and k! . The qubit transition frequency is u>o, while 
the constants g k represent the coupling between the k-th 
field mode and the qubit transition. As follows we refer 
to a TLS interacting with a cavity supplying a reservoir 
of field modes studied by Garraway [§] and adopted in 
Ref. [32]. 

Starting from an initial generic unentangled condition 
where the cavity is in the vacuum state |0)e, 

|*(0)) = (co|0)+ci(0)|l))®|0) B , 
the exact time evolution is described by the form 

oo 

|*(t)> =co\0)®\0) E + ci(t)\l)®\0) E +Y, b k(t)\Q)®\k)E, 

k=l 

where = a j-|0)_E for every k = 0,1,2, — The dy- 
namics is easily studied in the interaction picture, 

|¥(t))j = e< Hs+HB ^(t)) 

OO 

= cq\0) ® \0) E + Ci{t)\l) ® \0) E + B k (t)\0) ® \k) E , 
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ci(t), B k (t) = e lu "%(t), fc=l,2,. 



where i is the imaginary unit. The amplitude C\(t) is 
driven by the convoluted structure equation 



C 1 (t) = -(f*C 1 )(t), 



(1) 



where / is the two-point correlation function of the reser- 
voir of field modes, 

oo 
k=l 

For a continuous distribution of modes described by 
7](u>), the correlation function is expressed through the 
spectral density function J(uj), 

/>oo 

/(r)= / J(u)e-<"-^ T cbj, 



III. THE EXACT DYNAMICS 

We study the exact dynamics of the reduced density 
matrix of the qubit as given by Eq. @ for a reservoir 
described by the following spectral density function: 



J(w) 
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where O is the Heaviside step function, while A and a 
are arbitrary positive constants. Note that J (uS) ex- 
hibits a sub-ohmic behavior at low frequencies and has 
an absolute maximum, 3 3/4 A/ (2 a 3 / 2 ), at the frequency 

uim = + a/y/3. In contrast, the high frequency be- 
havior of J (u>) is an inverse power law like the usual 
Lorentzian one, though with a different power. The exact 
solution of Eq. Q corresponding to the spectral density 
([5]), can be shown to be 



G(t) 



^^Ti?(z,) z ie ^r(i/2,zfi), 



where R(z) is a rational function, 

(1-0 (a 1 ' 2 +z) (taW + z) 
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while the complex numbers z\,z%,Zs and £4, are the 
roots, distinct for every positive value of both A and 
a, of the polynomial Q(z), given by the following form: 



Q{z) 



f 2/aA + iaz 2 + + a 1/2 z 3 + z 4 



(8) 



For the sake of shortness, we do not report the analytical 
expressions of the roots, while the proof that the expres- 
sion ([6| is solution of Eq. Q is summarized in Appendix 
|A"1 The asymptotic expansion of the Incomplete Gamma 
functions [T3] in Eq. (|6j), identifies a time scale r and a 
decoherence factor V, depending on both the parameters 
A and a, 



t = max 



[\ Zl \ \ l = 1,2,3,4} 



(9) 



where J (lu) = r\ (lj) \g (w)\ and g (uj) is the frequency 
dependent coupling constant. 

The reduced density matrix, obtained by tracing over 
the degrees of freedom of the Hilbert space of the bosons, 
reads 



Pi,i(*) = l-A),o(*)=Pi,i(0) \G(t)f 



(2) 



pi,o(t)=Ph(t)=PiAO)e- Wot G(t). (3) 
The term G(t) fulfills the convolution equation: 

G(t) = -(/*G)(t), G(0) = 1. (4) 



V j_i 



(10) 



such that, for times f>r, a power law behavior emerges, 
described by the following asymptotic form: 



G(t) - -£>£- 3 / 2 , <^+c 



(11) 



We are finally equipped to give the exact time evolu- 
tion of the reduced density matrix through Eq. ^ and 
Eq. (§, 



Pi,i(0) 



^e z <*r(l/2,^) 



. (12) 
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obviously, p ,o(u = 1 - and = Pi,o(*)- 

Over long time scales, < 3> r, the dynamics is described 
by the asymptotic forms obtained through Eq. (Ill, 



Pi,i(*)~Pi,i(°)l 2? l 2 *~ 8 > i^+oo, (14) 



Pi,o(*)~-/»i,o(0)e- ,wo *2?r 3 / a , 



(15) 



Ultimately, the qubit collapses into the ground state. 

The dynamics of the qubit, driven by the special reser- 
voir described by the spectral density function J (uj), re- 
sults in a linear combination of Incomplete Gamma func- 
tions. In contrast to the exponential-like relaxation oc- 
curring for a Lorentzian spectral density function, in the 
present case the relaxation is Eulerian. The time scale r 
suggests the times for inverse power law relaxations with 
powers either 3/2 or 3, for either the coherence term or 
the populations, respectively. 




FIG. 1: The coherent term |pi,o(t)| of the reduced density 
matrix of a qubit interacting with a reservoir, described by ei- 
ther Jl (uj), both in strong (thin solid line) and weak coupling 
regime (dashed line), or J (uj) (thick solid line) spectral den- 
sity function, given by Eq. ( 13 1 and Eq. ( 19 1, respectively, for 



< t/r < 5.9. The values of the parameters are A — 0.8a 



5/2 

ujo — a/2, t ~ 0.974/a, tb = 1/a and tr — 1/(10 a) in strong 
coupling regime, tb = 1/(20 a) and tr = 10/(13 a) in weak 
coupling regime, the initial conditions read pi,i(0) = 1/2, 
Pi,o(0) = 1/5, and Co labels the cross point of the largest 
time coordinate between the thin and the thick solid line. 



IV. EXPONENTIAL VS EULERIAN 
RELAXATION 

It is now useful to compare the above results to the 
dynamics of qubit in a Garraway model [31 H2] j where 
the reservoir of field modes is described by a Lorentzian 
spectral density function 

^M = ~ 7 7A ! x2 - (16) 

2tt (uj - uj ) + A 2 

The corresponding bath correlation function is 

/*(*) = £e-*l«, (17) 

the positive constant A defines the spectral width of the 
coupling and is the inverse of the reservoir correlation 
time tb, while the positive constant 7 is the inverse of 
tr, the time scale on which the system changes. For 
details we refer to [3] . The time evolution of the reduced 
density matrix of the qubit reads 

Pi,i(t)=Pi,i(0)\G L (t)\ 2 , (18) 

Pi,o(t) = Pi,o(0)e- wt G L (t). (19) 
For A > 27, we have 

G L (t) = e-^ (cosh (| *) + 5 sinh (±t) ) , 
being d = \J A 2 — 27A, while, for A > 27, 



being d = y/2"/X — A 2 , with discrete zeros at times t n — 
2/d [n-K — arctan (d/\j ^ , for every n = 1, 2, 

We now compare the exponential-like decoherence pro- 
cess, given by a reservoir described by the Lorentzian 
spectral density function Jl (w), to the Eulerian relax- 
ation, resulting in an inverse power law over long time 
scales, corresponding to J (to), given by Eq. A 
long time scale inverse power law relaxation is slower 
than an exponential one, i.e., the decoherence process 
is strongly suppressed by adopting reservoirs described 
by J(oj), instead of Jl(lS). The qubit has "longer 
life". Numerical computations (Mathematica) for par- 
ticular values of the parameters are performed. For 
A = 0.8a 5 / 2 , wo = a/2, the following estimates emerge: 
zi a- 1 / 2 ~ -1.282 + i0.716, z 2 a' 1 ' 2 ~ -1.150 - 1 1.150, 
z^a- 1 /' 2 ~ 0.716-il.282, z 4 a~ 1 /2 ~ 0.717 + z 0.717, 
t ~ 0.974/a, \V\ ~ 0.112 a" 3 / 2 and the dynamics of the 
qubit can be evaluated starting, for example, from the 
initial conditions pi,o(0) = 1/5 and pi,i(0) = 1/2. 

For the damped Jaynes-Cumming model the absolute 
value of the coherence term shows an oscillating behav- 
ior enveloped in an exponential decay in strong coupling 
regime, 7 > A/2, and a linear combination of exponential 
decays in weak coupling regime, 7 < A/2, while, the Eu- 
lerian dynamics exhibits a power law decay proportional 
to t~ 3 l 2 over long time scales, i>T. A detailed behavior 
is shown in Fig. [T] for a short time scale, t < r, and in 
Fig. [2] for long time scales, f>r, respectively. 

In agreement with the exact analytical results ob- 
tained, the reservoir described by J (uj) gives an inverse 
power law relaxation, with power 3/2, over long time 
scales, t 3> r, given by Eq. |9]). If compared to the 
exponential-like decay corresponding to a reservoir de- 
scribed by J (lu), the decoherence process is slowed. 
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FIG. 2: The relaxation of coherent term |pi,o(t)| over long 
time scales, of the reduced density matrix of a qubit, interact- 
ing with a reservoir, described by either Jl (oj), both in strong 
(thin solid line) and weak coupling regime (dashed line), or 
J (u) (thick solid line) spectral density function, given by 
Eq. ([15} and Eq. jit)}, respectively, for 3.2 < t/r < 30. 
The values of the parameters are A — 0.8 a 5 ^ 2 , uio = a/2, 
r ~ 0.974/a, tb = 1/a and tr = 1/(10 a) in strong cou- 
pling regime, tb = 1/(20 a) and tr = 10/(13a) in weak cou- 
pling regime, the magnitude of the decoherence factor reads 
\V\ ~ 0.112a" 3/2 , the initial conditions are pi,i(0) = 1/2, 
Pi,o(0) = 1/5, and Co labels the cross point of the largest 
time coordinate between the thin and the thick solid line. 



V. CONCLUSIONS 

The exact dynamics of a TLS in the Jaynes-Cummings 
model, has been determined analytically by Garraway for 
Lorentzian type and a non-Lorentzian spectral density 
functions of the reservoir of bosons, in each case deco- 
herence results in an exponential-like relaxation. In line 
with the attempt to slow down the decoherence process, 
we design a reservoir piecewise similar to those usually 
adopted, i.e., a 1/2 power law behavior at low frequen- 
cies, similar to the sub-ohmic case, and a 3/2 inverse 
power law behavior at high frequencies, similar to the 
Lorentzian case, though with a different power. The spe- 
cial reservoir also exhibits a photonic band gap (PBG) 
edge O [Jjj] coinciding with the qubit transition fre- 
quency. 

The analytical description of the exact time evolution 
of the qubit is a combination of Eulerian functions of 
time and decoherence results in a 3/2 inverse power law 
relaxation over an evaluated long time scale. Ultimately, 
the system collapses into the ground state. If compared 
to the relaxation process corresponding to a Lorentzian 
type spectral density function, the process of decoherence 
is strongly suppressed. 

An environment implementing the specially designed 
reservoir of modes can in principle be realized with PBG 
media [Ml [17]. An anisotropic model providing a PBG 
close to the 3D photonic crystals [IB] , is discussed in Refs. 
[T5] and [TS]. The corresponding density of modes reads 



n (w) oc - oj e 6 ((oj - oj e ) /w e ) , 
where uj e is the band edge frequency. If the qubit tran- 



sition frequency coincides with the edge of the PBG, the 
low frequency behavior of the specially designed reservoir 
is recovered by assuming that the couplings vary slowly 
at low frequencies, g (oj) ~ g (ojo) for w ^ ojq. Physi- 
cally, the modes relevant for the dynamics are the reso- 
nant ones, which means that the time evolution mostly 
depends on the frequency behavior near the transition 
frequency of the system of interest [2U] . If a qubit can be 
placed in such a material and interact with such a reser- 
voir (in rotating wave approximation), the described de- 
coherence process can in principle appear. Of course, the 
accuracy of the model depends upon the high frequency 
behavior of the spectral density. 

Another theoretical model providing a structured PBG 
is the A^-period one dimensional lattice discussed in Ref. 
[21] , The density of frequency modes can reproduce a 
band gap by properly arranging the periodic sequence 
of unit lattice cells. The density of modes is evaluated 
analytically as a function of the complex transmission 
coefficients of each unit cell. 

Another potential realization of structured PBG are 
tunable ID PBG microcavities [35] [35] • Their fabrica- 
tion is achieved through advanced diffractive grating and 
photonic crystals technologies. 

The action of such structured PBG environments on a 
qubit could be a way of delaying the decoherence process 
with fundamental applications to Quantum Information 
Processing Technologies. In future we plan to search for 
specially engineered reservoirs inducing an even slower 
relaxation of qubit decoherence. 



Appendix A: the time evolution in details 

A detailed analysis of the equations driving the exact 
dynamics, is performed and the analytical solution of Eq. 
Q is evaluated in case the reservoir is described by the 
spectral density ([5}. 

The following class of non-negative, non-divergent and 
summable spectral density functions: 



J (uj) duj < oo, J (lj) = 6 ((uj — loq) /loq) A (oj — ojq) 



gives the Laplace transform 



G(u) = [u-tS(A)(-m)}- 1 , 



(Al) 



as a function of the Stieltjes transform S (A) (u), hold- 
ing true in case > and |arg{—iu}| < n. The 
fundamental uniform convergence of integrals involved 
is guaranteed by the constraints: L A (lj) du; < oo and 
> 0. We also notice that we adopt a positive range 
of continuous frequencies. The choice of the spectral den- 
sity function ([5} leads to the following Laplace transform: 



G(u) = 



a 1 / 2 da 1 ' 2 



i 1 ' 2 ) (a 1 ' 2 



,1/21 



(ttV2A + ia?l 2 u + (l + i) au 3 / 2 + o}l 2 u 2 ) 
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We are finally equipped to find G(t), solution of Eq. Q. 
The analysis of the discriminant shows that the roots of 
the polynomial Q{z), given by Eq. ([8]), are distinct for 
every value A > and a > 0, thus, the Inverse Laplace 
Transform of G(u), 

+-3/2 4 poo 

G(t) = —Y,R (zi) / re*'- T /Wdr, t > 0, 
v 71 " ;=1 Jo 

leads to Eq. ([6]). Also, the relations 

4 4 

Vi?(z ; )=0, lim G(i) = Zl = 1, 

obtained through the Initial value theorem, allow the 



continuation G(0) = 1. Thus, Eq. (|6| is the exact solu- 
tion of Eq. Q, where R(z) is given by Eq. Q. 

Finally, the asymptotic expansion of the incomplete 
Gamma function gives Eq. (11 1, describing the behavior 
of G(t) over long time scales, t ^> t, defined by Eq. (JoJ) , 
driving the dynamics over long time scale, as well. 
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